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F OREWORD 

This report, "High Speed Cylindrical Roller Bearing Analysis 
details information required to use the desigrt and analysis of 
computer program CYBEAN. The program described in this volume 
conforms to the 1978 version of CYBEAN prepared under U.S. govern 
ment Contract No. NAS3-20068 except that the program can now be 
executed at four levels of complexity. 

All efforts involved in the generation of this new version 
of the code were sponsored by the NASA-Lewis Research Center 
of Cleveland, Ohio, under the administration of the Bearing and 
Gear Analysis Section. The technical monitor was Mr. H. Coe. 

The work was performed under Contract'No. NAS3~22690 at SKF 
Industries, Inc., King of Prussia, Pennsylvania, during the 
period March 1981 through May 1981. 

Technical project leadership was executed by Mr. G. J. Dyba, 
with contribution from Mr. R. J. Kleckner. 
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HIGH SPEED CYLINDRICAL ROLLER BEARING ANALYSIS 
I. INTRODUCTION 

CYBEAN (Cylindrical BEaring ANalysis) has been created 
to detail radially loaded, aligned and misaligned cylindrical 
roller bearing performance under a variety of operating conditions. 
Emphasis has been placed on detailing the effects of high speed, 
preload and system thermal coupling. Roller tilt, skew, radial, 
circumferential and axial displacement as well as flange contact 
have been considered. Variable housing and flexible out-of-round 
outer ring geometries, and both steady state and time transient 
temperature calculations have been enabled. The complete range 
of elastohydrodynamic contact considerations, employing full 
and partial film conditions have been treated in the computation 
of raceway and flange contacts. The program properly directs 
roller-to-raceway contact load vectors automatically in those 
cases where roller and/or ring profiles have small radii of 
curvature. 

Volume 1 [7] of this report describes the models and associated 
mathematics used within CYBEAN. The user is referred to the 
material contained therein for formulation assumptions and algo¬ 
rithm detail. The material present in this. Volume II, is 
structured to guide the »ser in the practical and correct implement¬ 
ation of CYBEAN. Input and output architectures containing guide¬ 
lines for use and a sample execution are detailed in the matter 
whish follows. 

n. INPUT DATA 

CYBEAN requires the preparation of input data which in general 
describes the bearing geometry, properties of the materials used, 
and specifies the imposed operating conditions. With these inputs 
defined, optional solution procedure control may be selected. 
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The input architecture has been formulated to impose mini¬ 
mal initial demands on the user, Data is segregated into 
CATEGORIES which individually address specific characterization 
of the configuration addressed. Any single data i-et item 
required by the program falls into one of ten distinctly iden¬ 
tified subsets or categories (Table 1). Category "ROLLER", 
for example, contains all roller geometry data. Category "CAG" 
details cage information. Items detailing operating parameters 
such as load and speed are entered into input categories "0PER8" 
and “LOAD". 

All data required by the basic program are accepted in 
free NAMELIST format and default values are hardcoded to mini¬ 
mize the demands on user judgment. Special input data, required 
when the program options are used, is in 80 column card image 
foi’mat, 

In each category, the user has the freedom to specify 
all, part, or none of the data. If an item of data is omitted, 
a default value is assumed, A list of these default values 
is shown in Table 2. Failure to include basic data, e.g., ring 
radius, load, speed, etc, results in a diagnostic abort message, 

Dara comprising a category is specified in free format, 
with the restrictions that (!) column one of any card is not used 
and (2) all pieces of data in any category are separated by 
commas 1 . No specific sequence of data is required within each 
category. A minimum of two cards is needed to specify a complete 
set of data within a category. 

l 0n some computers the comma must follow the last significant. 
digit of an input variable. It is suggested that this restric¬ 
tion be observed to avoid the inconvenience caused by compiler 
peculiarities. 
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As an example, consider the category ’’IRING". Three items 
are needed. Using the nomenclature of Figure 1, they arei 

1. RIG - The groove radius of the inner ring 

2. FLGALI - The flange layback angle of a flanged 

inner ring on the left side. 

3. FLGAR1 - The flange layback angle of a flanged 

inner ring on the right side. 

An example of free format data is illustrated for the 
category ”IRING M in Figure 2, In this case the user wishes 
to describe the geometry of a flanged inner ring. Three cards 
are required. The first card contains a Dollar Sign ($)*in column 
2 followed by the category name "IRING". The second card is 
used to specify values of input data. Note that free format 
is used throughout and that all pieces of data are separated 
by a comma. The third card contains a Dollar Sign ($)*in column 
2 and the word END in columns 3 through 5 signifying the end 
of data for the category. Column 1 is never used in specifying 
data or category. 

Each data category is used to describe a particular aspect 
of either program use or the bearing configuration. Categories, 
in turn, must be arranged in the sequence noted in Table X before 
they can be used to transmit data to the program. 

The following paragraphs will list, in their proper order, 
all categories and the data they contain. In certain cases 
where, at the user’s option, categories can be omitted from 
the set, these options are made clear. Likewise, it is also 
clearly indicated when a category must be included in the set, 
regardless of execution options. If the user wishes to omit a 
category of data, he must still include the two cards: 

$CATEGORY NAME 
__ $END 

different computers may allow or require a different symbol. 
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CATEGORY 1 - SOLUTION CONTROL PARAMETERS 
CATEGORY NAME: SOLV_ 


CATEGORY DESCRIPTION { 

Computer program CYBEAN uses a Newton-Raphson 
iterative scheme to compute values for the governing 
equilibrium equation set, The user may wish to over¬ 
ride existing solution control parameters. Those, 
which are permitted as input by the user are con¬ 
tained within this category, 

DATA ITEMS; DEFAULT 

ITMAX - maximum number of iterations to be used in the 15 

Newton-Raphson iteration scheme, 


NPR - debug output print flag. Allows the user to see cal¬ 
culated results at intermediate steps of solution. NPR 
may be input with the following values; 


NPR * 0 No intermediate output printed 

NPR * 1 Divergence messages, intermediate 

equation residues (see Volume I) and 
roller-raceway loads are printed 

NPR 3 2 All output contained in NPR*1 plus 

the corrections of the variable values 
(see Volume I) as calculated in the 
linear equation solver 

NPR a 3 All output given for NPR=2 plus the 
solution algorithm used in this par¬ 
ticular execution 

NPR * 4 All output given for NPR=3 plus the 
matrix of partial derivatives 


skftechnol 
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DATA ITEMS DEFAULT 

CONVER - convergence criterion used to halt the iteration ^ 
procedure. Solution is said to be obtained when 

N 

( Z EQ.)/N < CONVER (1) 

i*l 1 

and 

It EU1EQ. | k - |EQ. | k ' 1 )/|EQ.| k }|< 2 (2) 

i* 1 x 

Here, k is the iteration index and EQ^ is the i-th equation 
residue. 

LEVEL - Solution control level: 1 

LEVEL =1: 

Inner ring and roller equilibrium is 
mutually satisfied through consideration 
of the elastic (Hertzian) contact loads 
at all roller-race contacts. Rollers are 
permitted to translate radially and rotate 
about their z-axis (roller tilt). The 
inner ring may translate in its radial 
plane. This execution level will permit 
economic calculation of bearing life. 

Lubricant and friction related loads are 
not considered; epicyclic roller speeds 
are assumed. 

LEVEL = 2: 

Mutual equilibrium of rollers and inner 
ring is sought as in Level 1, then used to 
evaluate the lubricant related forces (EI1D, 
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DATA ITEMS 

HD, drag, and cage pocket). This level o£ 
execution permits an economic estimation 
of bearing generated heat. Epicyclic roller 
speeds are assumed. 

LEVEL - 3!' 

Two solution algorithms are used to evaluate 
the bearing performance. The first seeks 
mutual equilibrium of inner ring and rollers 
by varying the inner ring position and roller 
radial, axial and tilt positions. The second 
algorithm seeks mutual equilibrium of the 
rollers and cage by varying the cage position 
and roller skew, tilt and rotational speeds. 
This level of execution provides an economic 
method of evaluating bearing performance when 
there is little or no rolling element skidding, 
and can be used to predict the existence of 
such skidding, but not its absolute magnitude. 
Note that at this level of execution, the 
inner ring reactive load may not exactly 
equilibrate the applied loading. 

LEVEL = 4: 

The two algorithms used in Level 3 are repeated 
until mutually satisfied. This assures equil¬ 
ibrium of the inner ring with applied loading. 
This execution level can be used to predict the 
presence and magnitude of rolling element 
skidding. 
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DATA ITEMS DEFAULT 

ITILT - Tilt analysis flag. T 

A .TRUE, value permits roller rotation 
about its z-axis (roller tilt) by solving the 
moment equilibrium equation. This flag is 
normally set to .TRUE.. However, in those 
cases where the roller tilt is expected to 
be very small, the user can specify ITILT * 
.FALSE, to save on execution time. 
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[CATEGORY 2 - PROGRAM LOGIC 
CATEGORY NAME: LOGIC _ 

CATEGORY DESCRIPTION : 

Within this category, the user is permitted to specify 
Values for logic used in a given program execution. User 
provided values dictate the program options. 

All variables in this category are "logical", and 
have either of two values, .TRUE, or .FALSE, eg., 

MPROP - .TRUE. 

In many cases, additional data will be required from 
the user as a consequence of selecting a specific program 
option. Descriptions of this extra input are found in the 
section "SPECIAL INPUT DATA", starting on page 23. 

DATA ITEMS: DEFAULT 


PLTRNG - Allows the user to obtain a line printer plot of 
the inner and outer ring profiles along their axial effec¬ 
tive length. 

PLTROL - Allows the user to obtain a line printer plot of 
the roller profiles along their axial effective length. 


ECHO - Allows the user to echo check the input data. This 
option invokes routines which print the data immediately 
after it has been read. 



F 


The following eight data items permit the user to 
invoke certain special program options. All items require 
the user to include additional data (see "SPECIAL INPUT DATA"). 
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DATA ITEMS DEFAULT 

COEF - Allows the User to input the influence coefficients p 
for the housing or other outer ring support structure. 

MPROP - Allows the user to input material properties for 
the rings, rollers, cage and housing. 

OVREND - Allows the user to input values of roller and p 

ring radii at ,the last 3 positions across the effective 
length. These values will overwrite values computed in 
the program. 

SYMY * A .TRUE, value indicates that the roller and rings 
are symmetric about their respective y-axis. If SYMY is 
.FALSE., the non-symmetric roller and race profiles must 
be read in. 

EVSLIC! - The program uses a slicing technique (see Volume I) 
to compute the roller-raceway loads. By default all slices 
are of equal width. In many cases it may be advantageous 
to specify that these slices are of unequal width. Unequal 
slice widths may be included as part of the input data by 
specifying EVSLIC = .FALSE.. 

FITS - Allows the user to use the clearance change portion p 
of the analysis. Default is no fit calculation, and a 
flexurally rigid outer ring is assumed. 

THERM - A .TRUE, value allows the user to use either the p 

steady state or time transient temperature calculating 
routines. 

00R - Used to denote an out-of-round outer ring. Operating 
diametral clearance for circular ring bearings may be cal- F 
culated with the program by specifying 00R = .FALSE, and 
FITS = .TRUE.. 
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CATEGORY 3 - ROLLER GEOMETRY DATA 
CATEGORY NAME; ROLLER 


CATEGORY DESCRIPTION 

Within this category the user must describe the 
geometry of the rolling elements within the bearing 
complement. This category is always included. 

(See Figures 3 and 4, all lengths are in mm.) 


DATA ITEMS 


DEFAULT 

ROLLD 

- Roller maximum diameter 

None 

RTL - 

Roller total length 

None 

RCR - 

Roller Crown Radius 

None 


SPHR - Roller end sphere radius on the right side of roller 33 ^ ^ 
(see Figure 3). 

SPHL - Roller end sphere radius on the left side of roller 381 mm 


RFL - Roller flat length None 

BLO -• Effective length of the roller outer ring contact 1 None 

ELI - Effective length of the roller inner ring contact 1 None 

XL, XR - The x-coordinate of the roller end sphere origin None 
(see Figure 3) 


l The effective contact length refers to the longest possible 
length which can be used to transmit load between the roller 
and raceway. Typically, this is the roller total length less 
corner radii. If, however, the raceway undercuts are excep¬ 
tionally large so that the track width is less than the roller 
effective length then the track width should be input. 
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DATA ITEMS DEFAULT 

BMPLAYO - Roller end-flange end play for a bearing having 
a flanged outer ring (see Figure 4). 

EMPLAYI - Roller end-flange end play for a bearing 

having a flanged inner ring (see Figure 4). 0 ' 

DIACL 1 - Bearing diametral clearance (see Figure 4), 

When using the FITS option (option 2) to calculate 0 

circular ring bearing operating clearance, DIACL is the 
cold, unmounted diametral clearance. 

KLUE « Roller geometry flag having the following four ^ 

possible values: 

KLUE ■ 1; The roller active profile is either fully 

flat at crowned with a flat. Symmetry about 
y is assumed. (KLUE“1 is applicable to the 
roller shown in Figure 3.) 

KLUE = 2; The roller active profile is fully crowned 
and symmetric about y. 

KLUE = 3; The roller active profile is symmetric about 
the y-axis but will be read in by the user 
(see "SPECIAL INPUT DATA"). 

KLUE = 4; The roller active profile is non-symmetric 
about the y-axis and will be read in by 
the user (see "SPECIAL INPUT DATA"). 

NUMROL - Total number of rollers in the complement (maximum ^ e 
permitted value = 50). 

■^Calculations with out-of-round components require DIACL to 
be calculated. See Option 2, page 25. 
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DATA ITEMS DEFAULT 

NS - Number of roller raceway slices used in the analysis. 5 
Note; If the user specifies symmetry about the roller 
y-axis then NS is the number of slices per symmetric half. 

If the user specifies no symmetry then NS is the total 
number of slices. 

PHI1 - Angular, location of the first element in the 
complement. In the nomenclature of Figure 5, the angle 
is measured CCW positive from the bearing y-axis. This 
input is in degrees. 
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CATEGORY 4: OUTER RING DESCRIPTION 
CATEGORY NAME: ORING 


CATEGORY DESCRIPTION : 

ORING is used to describe the geometry of the 
outer ring. Specification for lobing of the outer 
ring is made through the input of special data and 
will be discussed later. ORING is always included 
with input. Ring geometry is defined in Figure 1. 
Lengths and angles are to be specified in millimeters 
and degrees respectively. 


DATA ITEMS 

ROG - Groove radius of the outer race 

FLGALO - Flange angle of the flange located on the left 
side of the outer ring 

FLGARO - Flange angle of the flange located on the right 
side of the outer ring. If either FLGALO or FLGARO is 
left unspecified, then the outer ring is considered to 
be without flanges. 

DM - Pitch diameter 

KRING - Used to define the geometry class of rings. 

KRING=1; Both raceways have a flat profile (i.e. 

RIG=ROG= s O) and ring geometry is symmetric 
about the y-axis. 

KRING=2; Both raceways have a fully crowned profile 
(i.e. RIG^O and ROG^O) and their geometry 
is symmetric about the y-axis. 


DEFAULT 

o.(fiat) 
No flange 

No flange 

None 

1 
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KRING-3; Both raceway profiles will be read in as 
user specified input, symmetry is assumed 
(see "SPECIAL -INPUT DATA", Option 8, page 32 ). 

KRINGM; Both raceway profiles will be read in by user 
specified input, no symmetry is assumed (see 
"SPECIAL INPUT DATA", Option 10, page 34 ) . 
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CATEGORY 5 - INNER RING DESCRIPTION 
CATEGORY NAME; IRING__ 


CATEGORY DESCRIPTION ; 

Category IRING is used to describe the geometry 
of the inner ring. IRING is always included with 
input data. Ring geometry is defined in Figure 1, 

DATA ITEMS 

RIG - Groove radius of the inner ring 

FLGALI, FLGARI - Flange angles. See FLGALO and FIGARO 
in category ORING. 


DEFAULT 
0 .(flat) 
No flange 
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CATEGORY 6 - CAGE DESCRIPTION 
CATEGORY NAME; CAG 

CATEGORY DESCRIPTION ; 

The data items contained within this category 
are used to describe the geometry of the cage. This 
set of data must always be included. All lengths 
are in mm (see Figure 6). 

DATA ITEMS DEFAULT 

I RIDE - Cage type flag +1 

IRIDE * 1; the cage is inner ring land riding 

IRIDE *~1; the cage is outer ring land riding 

IRIDE * 0; the cage is rolling element riding 


RLDC - Rail-land diametral clearance None 
SRW - Single rail width None 
RLD - Rail-land diameter None 
CPCLR - Cage pocket radial clearance None 
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CATEGORY 7 - OPERATING CONDITIONS 
CATEGORY NAME; OPER8 _ 

CATEGORY DESCRIPTION ; 

Bearing operating conditions and operating tem¬ 
peratures are given in this category. Outer ring, 
inner ring and .flange temperatures are used to eval¬ 
uate the properties of the specified lubricant at 
these locations. Bulk temperature (BULKT) is used 
to evaluate the properties of the specified lubricant 
contained in the free space of the bearing cavity. 

This information is subsequently used in the calcu¬ 
lation of the viscous drag force acting upon the 
rolling elements. Data for this category must always 
be included. All temperatures are specified in degrees 
Celcius. 


DATA ITEMS (See Figure 7) 

SS - Shaft speed - RPM 

BULKT - Average temperature of lubricant in bearing cavity 

TRE - Rolling element temperature 

THSG - Housing temperature 

TSHFT - Shaft temperature 

TOR - Outer ring temperature 

TIR - Inner ring temperature 

TF1 through TF4 - Flange temperatures as shown in Figure 7 


DEFAULT 

None 

100 . 

100 . 

100 . 

100 . 

100 . 

100 . 

100 . 
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CATEGORY 8 - LUBRICATION DATA 
CATEGORY NAME; LUBE 


CATEGORY DESCRIPTION ; 

Within this category the user specifies lubri¬ 
cant properties and other data which relate directly 
Vo the lubricant or to the definition of friction 
related processes. 


DATA ITEMS 


NCODE - The user may specify particular lubricant proper¬ 
ties or simply select a value of 1 through 4 for NCODE, 

The latter selection obtains lubricant properties from 
a precoded table. Specific values 1 of NCODE and associated 
lubricant properties are shown in Table 3. The user may 
input lubricant properties not in Table 3 by specifying 
NCODE"0, 

ZTO, ZTI - Lubricant replenishment layer thickness 1 at the 
outer and inner rings, respectively, (mm) 

ZTFQ, ZTFI - Lubricant replenishment layer thickness at 
the outer and inner flanges, respectively, (mm) 


l At the present time the magnitudes of the inner and outer reple 
nishment layer thicknesses have not been correlated with flow 
rate, particular lubricants or bearing speed. The user is re¬ 
quired to establish proper values of the replenishment layer 
thickness. The following guidelines are suggested: 

1) To avoid starvation, replenishment layer thickness 
should be 1 or 2 times the EHD film thickness. 

2) Because of centrifugal force, intuition suggests the 
outer be thicker than the inner replenishment layer. 
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DATA ITEMS DEFAULT 

XGAV - Percent of lubricant occupying the bearing cavity? 5. 

0. < XCAV < 100. 

FRK - Lubricant friction coefficient, used in the Allen [1] 3 -07 

traction model. Typical values lie in the range 
0.05< FRK < 0,08. 

AKN - Computer program CYBEAN uses a model developed by 
Loewenthal [2] to compute EHD film thickness in point and 50. 

line contacts. The term AKN, the lubricant film thickness 
coefficient, appears in that equation. Typical values 
are 18. < AKN < 50. 

XMUCG - Coulomb friction coefficient used at the cage 
pocket-rolling element contact. If ZTO=ZTI=0, XMUCG is .0175 
applied. 

XMURC - Dry coefficient of friction at race contacts. If .0175 
ZTO=ZTI=0, XMURC is applied. 

XMUFL - Dry coefficient of friction at the flange contact. .0175 
If ZTFO=ZTFI=0, XMUFL is applied. 

The following data must be included if NCODE was specified 
as zero: 

VIS1 - Viscosity of lubricant (CENTISTOKES) at 100°F. None 

VIS2 - Viscosity of lubricant (CENTISTOKES) at 210°F. None 

RHQ60 - Density of lubricant (gm/cm 3 ) at X5 0 C. None 

G - Thermal coefficient of expansion (l/C°) None 

COND - Thermal conductivity (watts/M-Deg C) None 


2 As with replenishment layer thickness the amount of free lubri¬ 
cant should be correlated with the operating parameters. At 
this time such correlations do not exist. XCAV values of less 
than 5 percent are recommended. 

3 Numbers in brackets designate References listed in Section 6. 
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CATEGORY 9 - BEARING APPLIED LOAD 
CATEGORY NAME ; LOAD _ j 

CATEGORY DESCRIPTION : 

All bearing applied loads, either forces or 
misalignments are specified in this category. This 
data set must always be included (see Figure 5). 


DATA ITEMS DEFAULT 

FY - Radial load in Y direction (Newtons) 0* 

FZ - Radial, load in 1 direction (Newtons) 0* 

THETAZ « Misalignment about the z-axis (degrees) 0. 

THETAY - Misalignment about the y-axis (degrees) 0. 


) 
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CATEGORY 10 - SURFACE FINISH AND FATIGUE LIFE DATA 
CATEGORY NAME: LIFE ___ 

CATEGORY DESCRIPTION : 

Self-explanatory. This category must always 
be included. 


DATA ITEMS DEFAULT 

RMSROL - The RMS surface roughness of the roller (microns) .2032 
RMSIR - The RMS surface roughness of the inner ring (microns) .254 
RMSOR - The RMS surface roughness of the outer ring (microns) .254 
CIR - Life correction factor 1 for the inner ring !• 

COR - Life correction factor 1 for the outer ring 1* 


l The numbers input for CIR and COR are used to account for 
improved materials by multiplying the raceway fatigue lives as 
calculated by Lundberg-Palmgren methods. Typical life factor 
values for modern steels are in the range of 2 to 3. 

In the ASME Publication Life Adjustment Factors for Ball and 
Roller Bearings , the Material Factor D and the Material Process 
Factor E should be used multiplicatively as inputs for CIR and 
COR. 


The Lubricant Life Factor FL is calculated directly in the 
program and applies to both the inner and outer rings. 



m. SPECIAL INPUT DATA 
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The user activates extended program capabilities by invo¬ 
king up to a maximum o£ eleven options. Logic used to activate 
an option was detailed in the preceding material, and is summa¬ 
rized in Table 4. 

The user may employ as many options as necessary in a single 
program execution. The only restrictions are found in the input 
data sequence shown in Table 4 and in the use of specific input 
formats (Appendix A). 

All options require the user to specify additional infor¬ 
mation. This data follows immediately after the basic catego¬ 
rized data. 
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III.l OPTION 1; USER SPECIFIED MATERIAL PROPERTIES 

The user may specify the mate.rial properties of any bearing 
component. This is done by setting the logical variable 
MPROP * .TRUE. The appropriate properties are entered according 
to the card format shown in Figure Al. Unspecified variables 
(i.e., those set to zero or left blank) are assigned the following 
values: 

204083 N/mm 2 

. 3 

0.1124 x 10‘ 4 per °C 
7.806 gm/cm 3 


Modulus of Elasticity 
Poisson’s Ratio 

Coefficient of Thermal 
Expansion 

Density 


24 


AT81D049 


III.2 OPTION 2; FIT CALCULATIONS 


This option enables the user .to analyze cylindrical roller 
bearings manufactured with out-of-round outer rings. The user 
is given complete flexability in specifying the geometry by 
which preload is induced. The four most popular methods for 
generation of the latter with noncircular outer rings are shown 
in Figure 8. Input card format is specified in Figure A2. 

Variables which describe the shape of the noncircular com¬ 
ponents follow: 

MEAN RADIUS - This variable is defined in the nomenclature 
of Figure 9 and Figure 10 as 

„ _ d max + d min 

k MEAN 1 

It is to be noted that the mean radius must be specified 
for the raceway profile, outer ring outer surface profile 
and the housing profile (Figures 9 and 10). 

ECCENTRICITY RATIO - This variable is defined in the nomen¬ 
clature of Figure 9 and Figure 10 as 

e . 2D MAX . x 

(d max + d min) 

This variable is used to describe the magnitude of out-of- 
roundness manufacture into the raceway, ring outer surface 
and housing. 
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LODE ORIENTATION ANGLE - The lobe orientation angle, 0, 
is the angle measured clockwise positive from the bearing 
y-ajcis to the first lobe (Figures 9 and 10). Input is in 
degrees. 

NUMBER OF LOBES - Although most bearings which are manufac¬ 
tured out-of-round are made with 2 lobes, the user may 
input any number of lobes greater than or equal to zero. 

NOTE; Bearing diametral clearance under this option must 
be specified as follows: 

DIACL*{(MEAN RADIUS OF OUTER RACEWAY) - (RADIUS OF INNER 
RACEWAY SURFACE) - (ROLLER MAXIMUM DIAMETER) }x 2 

NOTE: The above calculations are with OOR * .TRUE.. Option 2 
can also be used for calculating the operating diametral clear¬ 
ance for circular ring bearings when OOR * .FALSE.. The input 
card format is specified on figure A2-A and the data cards 
would be placed in the same position as the OOR fit analysis 
cards. Definitions of items indicated for input card formats 
on A2-A are shown in Figure 11. 
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III.3 OPTION 3: USER SPECIFIED INFLUENCE COEFFICIENTS 


In the fit analysis described in Section III.2, it was 
assumed that the housing was rigid, therefore, all deformation 
was experienced by the outer ring. In some instances the user 
may wish to treat the outer ring support structure (the housing) 
as a deformable body. In this case, he must supply the ••influ¬ 
ence coefficients” 11 for the housing. The influence coefficients 
completely describe the deformation of the housing and are typi¬ 
cally obtained by using a finite element analysis. 

The user is required to input one coefficient per roller. 
Coefficients are evaluated at each roller location. Note that 
if the coefficients are not input, default is a rigid housing. 
The card image input format is shown in Figure A3. 


l The influence coefficients, Ci, are defined as the outward radial 
deformation of the housing at the i-th roller location due to 
a unit load at roller location 1. (min/N) 
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III.4 OPTION 4; USER SPECIFIED SLICE WIDTHS 

CYBEAN uses a slicing technique (Ref. 3 3 Ref. 7 of this 
report) to compute lubricant traction at the outer and inner 
ring roller contacts. Ordinarily, the slice widths are assumed 
equal, however, in some instances, it may be advantageous for 
the user to specify their individual and varying extent. It 
may, for example, be desired to obtain greater detail at the 
roller extremities. 

This option is invoked by specifying SYMY * .TRUE, and 
EVSLIC = .FALSE. Since symmetry is assumed, the slice widths 
need only be input for a symmetric half of the roller. The 
numbering scheme used is one where the first slice is encountered 
at the roller centerline, the last at the roller end (Figurel2), 
The card format shown in Figure A4 is used for input data. 
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III.5 OPTION 5; USER INPUT OF SYMMETRIC ROLLER GEOMETRY 

The specification of roller crown radius and flat length 
is sufficient to describe the geometry of the roller. If the 
user feels this to be insufficient, he may input values for 
the roller radii at specific locations along the roller 
effective length. The input sequence is shown in Figure 12. 

Note that the number of radii to be read in is equal to 
the number of slices plus one. Card formats are given in Figure 
A5. The specific example of logic used is SYMY * .TRUE, and 
KLUE - 3.. 
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III.6 OPTION 6: OVERWRITE CALCULATED fROLLER*) END RADII 
WITH USER SUPPLIED VALUES 

In cases where the roller effective length extends for 
nearly the complete roller length, it would be convenient for 
the user to specify the last 3 roller radii. In doing so, this 
would account for the. blend radius in the specification of the 
roller geometry. Again, note that symmetry is assumed, and 
the user need only specify the last 3 end radii on a symmetric 
half of the roller. 

The format of the input cards is shown in Figure A6, logic 
used is SYMY - .TRUE, and OVREND - .TRUE. 

This option also invokes Option 9. 
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111,7 OPTION 7i USER SPECIFIED. COMPLETELY VARIABLE 
ROLLER GEOMETRY 

The user may specify the complete detail of roller geometry 
as input by using the options: 

SYMY * .FALSE, and 
KLUE - 4 

When employing this option, the user must specify roller 
radii and slice widths across both the outer and inner ring 
effective lengths. Note that the number of slices input corres¬ 
ponds to the total number of roller raceway slices (Figure 13), 

Card formats for data input are shown in Figure A7. 
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iii,8 onion.. a.:. useb_sp_ecified symmetric ring geometry 

As with the roller geometry, the user is given the option 
of supplying the program with the symmetric ring geometry. The 
ring line of symmetry lies in the y-z plane of the bearing and 
divides the ring effective length into two equal parts. 

The data required when exercising this option is shown 
in Figure 14, input card formats in Figure A8, logic used is 
SYMY « .TRUE, and KRING - 3. 
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III.9 OPTION 91 OVERWRITE CALCULATED (RING) END RADII 
WITH USER SUPPLIED VALUES 

This option is similar to Option 6 except that the user 
specifies the last 3 ring radii instead of the last 3 roller 
radii. Note that when Option 6 is invoked, this option is also 
invoked, 

The format of the input data cards is shown in Figure A9, 
logic used is SYMY - .TRUE, and OVREND - .TRUE. 
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111,10 OPTION 10; USER SPECIFIED, COMPLETELY VARIABLE 
RING GEOMETRY 

The user may specify the ring geometry as input by using 
the options 

SYMY ■ .FALSE, and 
KRING - 4 

When using this option the user must specify all of the 
ring radii across both effective lengths. Note also, that the 
number of slices that are input corresponds to the total number 
of roller raceway slices (Figure 19- 

Card formats for data input are shown in Figure A10. 
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III.11 OPTION 11; TEMPERATURE CALCULATIONS 

CYBEAN may be used to compute either the time transient 
or steady state temperature distribution within a system defined 
by the bearing and its environment. Logic used requires 
THERM = .TRUE. 

The temperature portion of CYBEAN is designed to produce 
temperature maps for an axisymmetric mechanical system of any 
geometrical shape. The mechanical system is first approximated 
by an equivalent system which consists of a number of elements 
having simple geometric shape. Each element is then represented 
by a node point characterized by a mass, surface area, and having 
either a known or an unknown temperature. The environment sur¬ 
rounding the system is also represented by one or more nodes. 

With the node points properly selected, heat balance equations 
are formulated by the program for the nodes of unknown temper¬ 
ature. These equations become non-linear when there is radia¬ 
tion between two or more of the node points considered. 

The success of the approach depends largely on the realistic 
physical subdivision of the system. If the subdivision is too 
fine, there will be a large number of equations to be solved. 

If the subdivision is too crude, the results are likely to be 
inaccurate. 

The present thermal simulation is restricted to the treat¬ 
ment of axially symmetric physical systems. Bearing rings for 
example, fall into this category and can be represented by an 
element of uniform temperature. For a component or module which 
is not axially symmetric, the user must represent it with an 
equivalent axially symmetric element of approximately the same 
surface area and material volume. 

This section is based upon work presented in [4]. 
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With input data prepared as described in the following 
paragraphs, CYBEAN will solve the heat-balance equations for 
either the steady state or the time transient conditions and 
produce temperature maps for the physical system. 

INPUT DATA FOR TEMPERATURE CALCULATIONS 

Card formats for data input are listed in Figure All. 

Card 1 

Card 1 is a control card and contains input for both steady 
state and transient thermal analyses. It is not intended how¬ 
ever, that both analyses be executed with the same run. 

Item 1 : Highest Node Number (M). The temperature nodes 

must be numbered consecutively from one (1) to the highest node 
number. The highest node number must not exceed one hundred (100). 

Item 2 ; Number of Unknown Temperature Nodes (N). It is 
required that all nodes with unknown temperatures be assigned 
the lowest node numbers. The nodes which have known temperatures 
are assigned the highest numbers. 

Item 5 : Common Initial Temperature (TEMP)°C: The temper¬ 

ature solution iteration scheme requires a starting point, i.e., 
guesses of the equilibrium temperatures. Card 2 allows the user 
to input guesses of individual node temperatures, however, when 
a node is not given a specific initial temperature, the tempera¬ 
ture specified as Item 3 of Card 1 is assigned. 
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Item 4 : Punch Flag (IPUNCH): If thePunch Flag is not 

zero (0) or blank, the system equilibrium temperatures along 
with the respective node numbers will be punched according to 
the format of Card 2. This option is useful if, for instance, 
the user makes a steady state run with lubrication, and then 
wishes to use the resultant temperature as the initiation point 
for a transient dry friction run in order to assess the conse¬ 
quence of lubricant flow termination. 

Item 5 : '’Output Flag" (IUB). If the "Output Flag" is 
not zero the bearing program output and a temperature map will 
be printed after each call to the bearing solution scheme. This 
printout will allow the user to observe the flow of the solution 
and to note the interactive effects of system temperatures and 
bearing heat generation rates. Since the temperature solution 
is not mathematically coupled to the bearing solution the pos¬ 
sibility exists that the solution may diverge or oscillate. 

In such a case, study of the intermediate output produced by 
the "Output Flag" option may provide the user with better initial 
temperature guesses that will effect a steady state solution. 

Two levels of bearing output are permitted. If IUB is 1, the 
rolling element output is not required. If IUB is 2, full 
bearing output is obtained. 

Item 6 ; "Maximum Number of Calls to the Bearing Program" 
(IT1). IT1 is the limit on the number of Thermal-Bearing iter¬ 
ations, i.e., the external temperature equilibrium calculation. 
The user must input a non-zero integer such as 5 or 10 in order 
for CYBEAN to iterate to an equilibrium condition. If IT1 is 
left blank or set to zero (0) or 1, bearing performance will 
be based on the initially guessed temperatures of the system. 
Temperatures printed will be based on the bearing generated heats. 


37 




AT81DQ49 


It is unlikely that an acceptable equilibrium condition will 
be achieved. However, the temperatures which result may provide 
better initial guesses for a subsequent run than those specified 
by the user. 

IT1 also serves as a limit on the transient temperature 
solution scheme, by limiting the number of times the bearing 
solution scheme is called. Each call to the bearing scheme will 
input a new set of bearing heats to the transient temperature 
scheme until a steady state condition is approached or until 
the transient solution time-up limit is reached. 

Item 7 ; "Absolute Accuracy of Temperatures for the Exter¬ 
nal Thermal Solution" (EPI). In the steady state thermal solu¬ 
tion scheme, each calculation of system temperatures occurs 
after a call to the bearing scheme which produces bearing gene¬ 
rated heats. After the system temperatures have been calculated 
for each iteration, using the internal temperature solution 
scheme, each node temperature is checked against the nodal tem¬ 
perature at the previous iteration. 

If {t^^ “ (N-l) i^ <EP1 ^ or a11 nocies i then equilibrium 

has been achieved and the iteration process stops. 

Item 8 : "Iteration Limit for the Internal Thermal 
Solution" (IT2). After each call to the bearing program, the 
internal temperature iteration scheme is used to determine the 
steady state equilibrium temperatures based on the calculated 
set of bearing heat generation rates. If IT2 is left blank 
or set to zero (0), the number of internal iterations is limited 
to twenty (20). 
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Item 9 : "Accuracy for Internal Thermal Solution" (EP2). 

The use of EP2 is explained in Volume I. If EP2 is left blank 
or set to zero (0), a default value of 0.001 is used. 

Item 10 : "Starting Time" (START) is a time at which the 
transient solution begins, T g ; usually set to zero (0). 

I tern 11 : "Stopping Time" (STOP) is the time in seconds 
at which the transient solution terminates, T£. The transient 
solution will generate a history of the system performance which 
will encompass a total elapsed time of 


(Tf - T ) seconds ' 

Item 12 : "Calculation Time Step" (STEPIN). The transient 
internal solution scheme solves the system of equations (see 
Volume I): 


t k+l 



T = STEPIN 

The user may specify STEPIN. If left blank or set to zero 
(0), CYBEAN calculates an appropriate value for STEPIN using 
the procedure described in [4]. 

Item 13 : "Time Interval Between Printed Temperature Maps" 
(TTIME) seconds. The user must specify the length of time which 
will elapse between each printing of the temperature map. The 
interval will always be at least as large as the "calculation 
timestep" (STEPIN). 

Item 14 : "Time Interval Between Calls of the Bearing 
Program" (BTIME). BTIME will always have a value larger than 
or equal to (STEPIN) even if the user inadvertently inputs a 
shorter interval. Computational time savings result if BTIME 
is greater than STEPIN, however, accuracy might be lost. 
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Card 2 

In the steady state analysis this card is used to input 
initial guesses of individual nodal temperatures for unknown 
nodes as well as the constant temperatures for known nodes, such 
as ambient air and/or an oil sump. 

In the transient analysis. Card 2 is used to input the 
nodal temperatures of all nodes at time = T g , i.e., at the ini¬ 
tiation of the transient solution. 

Card 5 

With this card, node numbers are assigned to the components 
of the bearing. With this information the proper system temper¬ 
atures are carried into the respective bearing analysis. The 
inner race and inner ring node numbers may or may not be the 
same at the user’s discretion. Similarly, the Outer race and 
outer ring node numbers may or may not be the same. The flange 
numbering scheme is shown in Figure 1. 

Card 4 

The beaming analysis accounts for frictional heat generated 
at five locations in the bearing, i.e., the inner race, the 
outer race, between the cage rail and ring land, the bulk lubri¬ 
cant due to drag and at the flanges. The heat generated at the 
cage-rolling element contact is added to the bulk lubricant. 

This card allows the heat generated to be distributed equally 
to two nodes. For instance, the heat generated at the inner 
race-rolling element contact should be distributed half to the 
rolling element and half to the inner race. The heat developed 
between the cage and inner ring land may be distributed half 
to the inner ring and half to the cage if a cage node has been 
defined otherwise, half to the bulk lubricant. 
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Card 5 

This card specifies the node numbers and the heat generation 
rate at those nodes. This card is used to specify where heat 
is generated at a constant rate such as at rubbing seals or 
gear contacts. 

Card 6 

This card type is used to input the numerical values of 
the various heat transfer coefficients which appear in the equa¬ 
tions for heat transfer by conductivity , free convection , 
forced convection , radiatio n and fluid flow . Up to ten coeffi¬ 
cients of each type may be used. Separate values of each type 
of coefficient are assigned an index number via card 6 and in 
describing heat flow paths (Card 7 below) it is necessary only 
to list the index number by which heat transfers between node 
pairs. 

Indices 1-10 are reserved for the conduction coefficient 
X, 11-20 for the free convection parameters, 21-30 for forced 
convection , 31-40 for emissivity and 41-50 for fluid flow (pro¬ 
duct of specific heat, density and volume flow rate). 

As an example, for heat transfer by conduction with coeffi¬ 
cient X of 53.7 watts/M°C one could prepare a card 6 with the 
digit 1 punched in column 10 and the value 53.7 punched in the 
field corresponding to card columns 11-20. If a conduction 
coefficient of 46.7 were applicable for certain other nodes in 
the system one could punch an additional card assigning index 
No. 2 to the value X = 46.7 by punching a "2" in card column 
10 and 46.7 anywhere within card columns 11-20. 
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Rather than inputting constant forced convection coeffi¬ 
cients, optionally, these coefficients can be calculated by the 
program in one of three ways. If the calculation option is 
exercised a pair of cards is used in place of a single card 
containing a fixed value of a. The contents of the pair of 
cards depends upon which of the three optional methods are used. 

Option 1) a is independent of temperature but is calcu¬ 
lated as a function of the Nusselt number which 
in turn is a function of the Reynolds number 
R g , the Prandtl number P r as follows, (cf.[5]): 

«* 

Nu = KR* ? h r 

where X o; q is the lubricant conductivity, L is 
a characteristic length (with the units of meters) 
and K, a and b are constants. 

Option 2) a is a function only of fluid dynamic viscosity 
and viscosity is temperature dependent. 

ct = c 

Option 3) a is a function of the Nusselt, Reynolds and 
Prandtl numbers and viscosity is temperature 
dependent. 

Appendix B has been included to aid the user in data pre¬ 
paration and calculation of heat transfer coefficients. 
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Card 7 

This card defines theheat flow paths between pairs of nodes. 
Every node must be cr nnected to at least one other node, i.e,, 
two or more independent node systems may not be solved with a 
single program execution. 

The calculation of heat transfer areas is based on lengths, 
and L^ input using card 7. Additionally, the type of surface 
for which the area is being calcuated is indicated by the sign 
assigned to the heat transfer coefficient index. If the surface 
is cylindrical or circular the index should be positive, if the 
surface is rectangular the index should be input as a negative 
integer. 

In the case of radiation between concentric axially symme- 
trie bodies, is the radius of the larger body. For radiation 
between two parallel flat surfaces or for conduction between 
nodes, Lj is the distance between them. 

Fluid flow heat transfer accounts for the energy which the 
fluid transports across a node boundary. Along a fluid node 
at which convection is taking place, the temperature varies, 

The nodal temperature which is output is the average of the 
fluid temperature at the output and input boundaries. If the 
emerging temperature of the fluid is of interest, it is neces¬ 
sary to have a fluid node at the fluid outlet. At this auxili¬ 
ary node only fluid flow heat transfer occurs and the fluid 
temperature would be constant throughout the node. Thus the 
true fluid outlet temperature will be obtained. 

Conduction of heat through a bearing is controlled by index 
51. The actual heat transfer coefficient which contains a con¬ 
ductivity, area and a path length term is calculated in the 
bearing portion of the program. The term is based upon an aver¬ 
age outer race and inner race rolling element contact. 
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Card 8 

This card inputs data required to calculate the heat capa¬ 
city of each node in the system.' This card type is required 
only for a transient analysis. 
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IV. OUTPUT data 

CYBEAN output data is structured such that .all user supplied 
information, e.g. bearing geometry and operating conditions, are 

printed first. Default values assigned for unspecified input 
parameters are also displayed. Computed data resulting from 
the bearing analysis are then presented. Two examples are 
il uiUt ad representing solution levels 2 and 3 to illustrate typical 
program output. The applicable cylindrical roller bearing 
geometry and operating conditions are specified in Table 5. 

The following program options were used in the sample 
executions: 


PLTRNG 

= 

.FALSE 

PLTROL 

= 

.FALSE, 

MPROP 

= 

.TRUE. 

THERM 

= 

.FALSE, 

FITS 

= 

.FALSE, 
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IV.1 SAMPLE PROBLEM OUTPUT DESCRIPTION 

A* Output at Solution Level » 2 

The output generated at solution level 2 is presented in 
Appendix D. Output contained on the first page informs the 
user of the program version, and references the latest user's 
manual. User invoked options for this execution are also noted 
on this page. 

Pages D1 and D3 present an organized list of the basic 
categorized input data. Roller and ring geometry are listed 
on Page D4. 

Cage geometry, material properties and friction data are 
listed on Page D5 along with bearing component temperatures. 

The bearing fatigue life, as well as individual L 1Q fatigue 
lives of the outer and inner rings are presented on Page D6. 

The bearing life represents the statistical combination of the 
two raceway lives. The raceway lives in turn reflect the 
combined effects of the user input material factors and the 
computed lubricant film thickness factors. The film thickness 
to surface roughness ratio is used in the calculation of the 
lubricant life factor. Detailed information for this 
calculation is given in [6]. 

Pages D7 and D8 list the roller raceway contact loads at 
the outer and inner ring contacts as defined in the R coordinate 
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frame, illustrated in Figure 16. These forces include both 
elastic and lubricant traction effects. Hollers are numbered 
in ascending order beginning with the roller lying on (or nearest) 
the bearing y-axis and proceeding counter clockwise. On Page D9, 
the inner ring applied forces, moments and displacements, which 
constitute the system loading information supplied by the user 
are displayed. The calculated "inner ring reactive forces and 
moments" are the resultant of the vector sum of all roller-inner 
ring contacts. Load values are included with the output to 
assess the degree of convergence at the lower solution levels. 

In this LEVEL 2 execution, the user specified applied load was 
4448N in the radial Y direction. The program arrived at a solution 
of: 

Y * 4450 N 
Z = 1.69 N 

The theoretically correct solution requires Z “ 0, and the 1.69N 
residual represents only 0.04% deviation from the applied load. 

The lubricant data shown on the same page is self-explanatory. 
Output contained on Page DIO is for the most part self-evident. 

The roller orbital and rotational speeds printed are epicyclic 
since at level 2 they are not solved for. The pocket loads are 
estimated from rolling element Z equilibrium force equation 
residues. 

The Hertzian contact stresses provided on Page Dll represent 
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maximum values for the line (roller-raceway) and point (roller 
end-flange) contacts. At this solution level flange-roller 
contact is not analyzed. The lubricant film thicknesses on 
Page D12 represent minimum values for the line (roller-raceway) 
contacts. On Page D13, roller skew and tilt are presented to 
the user as measured in two distinctly different reference 
frames. [Note: At LEVEL 2, roller skew is not calculated.) 

Calculated roller skew and tilt is presented to the user as 
measured in two distinctly different reference frames. "Absolute" 
refers to the rotation the roller experiences relative to its 
initial position. "Relative" refers to the rotation the roller 
experiences relative to the inner ring position. The following 
example illustrates the two conventions. 

Consider a bearing whose four rollers are "frozen" in their 
position of zero absolute skew and tilt. Assume rollers to be 
located at 90° intervals, two being on the Y-axis and two on 
the Z-axis. The inner ring is now rotated about the Z-axis. 

With the ring in its final position, the rollers still have zero 
"absolute" skew, however, the rolleu which lie on the axis of 
rotation appear skewed when viewed from the "relative" reference 
frame of the inner ring. 

Heat generation rates are displayed on Page D14 for 
bearing components. These data are self-explanatory. 
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B. Output at Solution Level ~ 3 

The output generated at this solution level is presented 
in Appendix E. Data listed on Pages El through E8 are displayed 
in the same manner as for solution level 2, described in 
Section IV.1 A. 

On Page E9, in addition to applied and calculated inner ring 
reactive forces, misalignment of the inner ring was considered. 

At LEVEL 3 flange-roller contact can be analyzed. Flange 
induced roller loads are listed on Pages E10-E11. The numbering 
scheme is identical to the one used for roller-raceway loads. 

These loads are defined in the E coordinate frame of Figure 16. 
Page E12 displays the sliding spefcs*; magnitude at the roller end- 
flange contact. The lubricant data is also represented on this 
page. At level three the roller speeds are computed using the 
roller equilibrium equations considering friction forces and 
moments. 

On Page E13, the calculated roller speeds are provided. 
Epicyclic speeds are printed for the user's reference. 

The cage pocket normal force listed is that force experienced 
by the roller due to interaction in the z direction with the 
cage web. A negative sign indicates that the roller is pushing 
the cage. All pocket loads are computed such that the rolling 
elements and cage are in equilibrium. 
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Maximum Hertzian stresses are displayed/ on Page E14 for 
the roller-raceway and roller-end flange contacts. The lubricant 
film thicknesses printed on Page BIS represent minimum values 
calculated for both raceways and the flange-roller contact. 

Page E16 lists the roller tilt and skew angles. Roller 
skew is calculated at this level (LEVEL 3). Roller tilt was 
not considered in this particular analysis (ITILT was set to 
.FALSE.). Thus, the tilt values presented on this page are not 
calculated, rather they are estimated by the program. 

Page E17 provides the heat generated at the rolling element- 
flange contacts in addition to the other bearing heat generation 
rates. 
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V. PROGRAM LIMITATION AND SPECIAL CASES 

CYBEAN is a design tool. As with any tool, successful 
use requires awareness of intended applicability and inherent 
limitations. 

A. LIMITATIONS 

The user must conform to the following geometric and oper¬ 
ating restrictions: 

1) The bearing complement may contain no more than 
fifty (50) rollers. However, this version of 
CYBEAN has been dimensioned to accommodate a 

203 x 203 array of partials [COMMON/BIGC/C£.41209)), 
This corresponds to only 28 rollers. The rest of 
the program is structured to accept up to 50 rollers. 
To execute CYBEAN with more than 28 rollers, the 
C array needs to be expanded. The equation to 
use in computing the array size is: 

DIMENSION = [7 X NUMROL + 7] 2 
NUMROL - NO. OF ROLLERS 

The /BIGC/ common appears in the following subroutines 

PARDA, PARDER, SOLV14, SIMQA, WRITC, 

MINT, CINT, RINT, RINFIT 

2) Flanges may be specified on the outer or inner 
ring, but not on both simultaneously. 

3) Given a cylindrical roller bearing operating 
with specified misalignment, and/or no geometric 
symmetry, one ring (either the inner or outer) 
must be flanged. 
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4) This edition of CYBEAN does not accept externally 
imposed axial loads. 

5) Extremely light radial loading*, wherein a single 
roller interacts with the inner ring, will cause 
error termination. 

6) Use of the fit option (i.e. FITS = .TRUE.), requires 
the mean radius of the outer ring outer surface to 

be less than the mean radius of the housing (for 00R = .TRUE. 


B. SPECIAL APPLICATIONS 

Some special applications of CYBEAN are: 

1) The user may approximate the bearing heat generation 
rate by specifying ITMAX=1. Heat generation rate 
computations under this option are based upon a single 
iteration of the initial guess independent variable 
values. 

The initial guess variable values are obtained by 
solving the governing equation set for the equilibrium 
of elastic forces. 

2) Through the use of the symmetric ring and roller geo¬ 
metry program options, the user can use the cylindri¬ 
cal roller bearing program to analyze a single row 


*A current estimate used to determine the minimal radial load is 


Here, C is the basic (AFBMA) dynamic capacity and P . is the 
minimum radial load. min 
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spherical roller bearing. When doing so, ring misa¬ 
lignment must be set to zero, 

3) The user may, when making several steady state temper¬ 
ature program executions, use the card punch option 
(IPUNCHj/0) to obtain the temperatures in 80 column 
card format. These provide an economic initial guess 
(see: Temperature Calculations, Card 2) for subsequent 
runs. 

4) The program is capable of properly directing roller- 
to-raceway contact load vectors automatically, in those 
cases where roller and/or ring profiles have small 
radii of curvature. 
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FIGURE 1 : USER INPUT RING DATA. 
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FIGURE 4 : USER INPUT BEARING CLEARANCES 










FIGURE 6 : USER INPUT CAGE GEOMETRY. 
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FIGURE 7: TEMPERATURE NODE IDENTIFICATION SCHEME. 
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FIGURE 8: USER SPECIFIED INFORMATION REQUIRED FOR FOUR MOST 
POPULAR MODES OF INDUCING ROLLER PRELOAD. 
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FIGURE 11: INPUT GEOMETRY DEFINITIONS FOR CIRCULAR RING FIT ANALYSIS 
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FIGURE 12: OPTIONAL SYMMETRIC ROLLER GEOMETRY INPUT DATA 








Line of 
Symmetry 



FIGURE 14: OPTIONAL SYMMETRIC RING GEOMETRY INPUT DATA. 



Effective 
Leu gth 



USER'S NOTE: 

THE NUMBER OF RADII INPUT IS EQUAL TO THE NUMBER OF SLICES PLUS ONE. 

FIGURE 15: OPTIONAL USER INPUT RING GEOMETRY. 
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FIGURE 16: ROLLER COORDINATE 
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FLANGED INNER RING 
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ORGANIZATION 

OF INPUT DATA CATEGORIES 


CATEGORY 

NUMBER 

CATEGORY 

NAME 

CATEGORY DESCRIPTION 

NECESSARY 

X 

SOLV 

Solution Control Parameters 

No 

2 

LOGIC 

Program Control Logic 

No 

3 

ROLLER 

Roller Geometry Data 

Yes 

4 

ORING 

Outer Ring Description 

Yes 

5 

I RING 

Inner Ring, Description 

Yes 

6 

CAG 

Cage Description 

Yes 

7 

OPER8 

Operating Conditions 

Yes 

8 

LUBE 

Lubrication Data 

Yes 

9 

LOAD 

Bearing Applied Loads 

Yes 

10 

LIFE 

Fatigue Life Data 

Yes 
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PROPERTIES OP FOUR LUBRICANTS 
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ROLLER BEARING SPECIFICATIONS 


Inner Race 


Bore Dia. 

mm (in) 

118 

(4.6457) 

Raceway Dia, 

mm (in) 

131.66 

(5.1834) 

Flange Dia. 

mm (in) 

137.47 

(5.4122) 

Width 

mm (in) 

26.92 

(1.060) 

Groove Width 

mm (in) 

14.59 

(0.5746) 

Flange Angle 


.6 i 

deg. 


Outer Race 


Outer Dia. 

mm 

(in) 

164.65 

(6.4760) 

Raceway Dia. 

mm 

(in) 

157.08 

(6.1842) 

Width 

mm 

(in) 

23.9 

(0.942) 


Rollers 


Diameter 

mm (in) 

12.65 

(0-4979) 

Length - overall 

mm (in) 

14.56 

(0. 5733) 

- effective 

mm (in) 

13.04 

(0.5133) 

- flat 

mm (in) 

8.40 

( 0. 3307) 

Crown Radius 

mm (in) 

622.3 

(24.5 ) 

End Radius 

mm--(in) 

381.0 

(15.0 ) 


Number 28 


Ca ge 

Land Dia. mm (in) 137.95 (5.4312) 
Axial Pocket Clea-rance mm (in) .020 (0.0008) 
Tangential Pocket Clearance mm (in) .221 (0.0087) 
Single Rail Width mm (in) 4.6 (0.18 ) 


Operating Conditions 

Shaft Speed 
Bearing Radial Load 
Oil Inlet Temperature 

Misalignment of Races 

Lubricant 


20,000 rpm 

4450 N (1000 lb) 

366.5K (200°F) 

0 deg. and 0.0833 deg. 
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FIGURE Al: USER INPUT MATERIAL PROPERTIES CARD FORMAT 






















FIGURE A2: FIT CALCULATION USER INPUT CARD FORMAT (PAGE I of 2) 
CARD 1 























USER input CARD FORMAT 












USER INPUT SLICK WIDTH CARD FORMATS 
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RDS AS NEEDED, FOLLOWED 3Y A BLANK CARD 






CARDS/COEFFICIENT, AS MANY AS NEEDED, FOLLOWED BY A BLANK CARD, 



FIGURE All (CONTINUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION. 










CARD 7 - HEAT FLOW PATHS; USE AS MANY CARDS AS NEEDED, FOLLOWED BY A BLANK 



FIGURE All (CONTINUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION- 








CARD 8 - NODAL HEAT CAPACITIES (USE ONLY FOR TRANSIENT CALCULATIONS) 



FIGURE All (CONTINUED): INPUT. DATA CARD FORMATS FOR TEMPERATURE CALCULATION. 
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B.l BASIC EQUATIONS* 

B,1.1 Heat Conduction 

The rate of heat flow q c ^ j(W) that is conducted from node 
i to node j may be expressed by, 


q 


ci, j 


t*i ■ V 

ij 


t. and t. are the temperatures at i and j, respectively, A 1 . the 

i j 2 x f •> 

area normal to the heat flow, (m ) the distance (m) and 
the thermal conductivity between i and j, (W/m°C). 


Assuming that the structure between point i and j is com¬ 
posed of different materials, an equivalent heat conductivity may 
be calculated as follows: 



The calculation of the areas will be discussed in Section B.l.5. 
*This Appendix is based on the material present in Reference 4. 
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B, 1,2 Convection 

The rate of heat flow that is transferred between a solid 
structure and air by free convection may be expressed by 


q vi,j 

where 


A i,i 




i,25 


SIGN (t i 


SIGN - 


( 1, if (t^ 

(-1, if (tj 


t j ) > 0 

V < 0 


V 


in which 




'2.5 • 10 ^ W/m 2 - (degC) 1,25 for hot surfaces facing 
upward and cold surfaces facing downward 

« - 1.4 • 10" 2 W/m 2 - (degC) 1,25 for hot surfaces facing 
downward and cold surfaces facing upward 

yL,8 • 10“ 2 W/m 2 - (degC) 1,25 for vertical surfaces 


For other special conditions, a^. must be estimated by referring 
to heat transfer literature. 

The rate of heat flow that is transferred between a solid 
structure and a fluid by forced convection may be expressed by 




°i,j A i,j Ct i • V 


in which ct-j is the convective heat transfer coefficient. 
Now, with a - , introduce the Nusselt number 

N - 
U A 

the Reynolds number 

R - UL 

R e " “ 

and the Prandtl number 

p = 
r A 
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where 

L is a characteristic length which is equal to the diameter in 
the case of a cylindrical surface and is equal to the plate 
length in case of a flat surface (m), 

U is a characteristic velocity which is equal to the difference 
between the fluid velocity at some distance from the surface 
and the surface velocity (m/sec). 

X is the fluid thermal conductivity (W/M°C) 
v is the fluid kinematic viscosity (M /sec) 
p is the fluid density (Kg/m’*) 

C n is the fluid specific heat (J/kg°C) 

For given values of R e and P r the Nusselt number N u and 
thus the heat transfer coefficient may be estimated from one of 
the following expressions: 

Laminar flow along a flat plate: R g < 2300 

n u ■ 0.323 /k; • 3 /r^ 

Laminar flow of a liquid in a pipe: 

N u ■ 1.36 VR e • 

where D is the pipe diameter and L the pipe length 

Turbulent flow of a liquid in a pipe: 

N u a 0.027 • Rg' 8 * 3 /p7 

Gas flow inside and outside a tube: 

N a 0 . 3 R°' 5 7 
u e 

Liquid flow outside a tube: 
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Forced convection from the outer surface of a rotating 
shaft 

N u - 0,11 [ 0.5 R e 2 • P r ] 0,35 

where the Reynolds number R e is developed by the shaft 

r . P 2 
e v 

in which to is the angular velocity (rad/sec) 

D is the shaft diameter (m) 

The average coefficient of forced convection to the lubri¬ 
cating oil within a rolling contact bearing may be approximated 
by, 

a - 0.0986 [1 + H L S 22 . JM. ] ) l »x(P ) 1/3 

using + for outer ring rotation 
- for inner ring rotation 

in which N is the bearing operating speed (rpm) 

D is the diameter of the rolling elements (mm) 
d m is the bearing pitch diameter (mm) 

$ is the bearing contact angle; zero for cylindrical 
roller bearings (degrees) 

B.1.3 Fluid Flow 

The rate of heat flow that is transferred from fluid node 
i to fluid node j by fluid flow is 

“fi.i ■ ^5 c p ct i * V 

is the volume rate of flow from i to j. It mus.t be observed 
that the continuity of mass requires the following equation to 
be satisfied 



provided the fluid density is constant. The summation should be 
extended over all nodes i within the fluid which have heat 
exchange with node j by fluid flow. 
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B.X.4 H eat Radiation 

The rate of heat flow that is radiated to node j from 
node i is expressed by 

«RiJ ' «t. ♦ 273D 4 - (t. * 2?i) 4 > 

where 

T. - t. + 273.16 
3 3 , 

T i * ti + 273.16 

and the value of the coefficient <$. . depends on the geometry and 

i > 3 

the emissivity.or the absorptivity of the bodies involved. 

Fo, ,f radiation between large, parallel and adjacent sur¬ 
faces of equal area, A. . and emissivity, e. *, 6. . is obtained 
" i>3 i >j 1 >3 

from the equation 


i,3 


,3 aA i,3 


where a, the Stefan-Boltzmann constant, is 
a = 5.76 • 10‘ 8 W/m 2 /CdegK) 4 


For radiation between concentric spheres and coaxial 
cylinders of equal emissivity, e^ j is given by the 

equation 


4,i aA i,i 
5T~ 




where a is as above A^ j is the area of the enclosed body and 

A*. . is the area of the surrounding body, i.e., A. . < A*. .. 

1,3 i,J X ,J 

Expressions for 6. . that are valid for more complicated 

1,3 

geometries or for different emissivities may be found in the 
heat transfer literature. 
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B.1,5 Calculation of Areas 

In the case of heat transfer in the axial direction A^ j 
is given by the equation (Fig* B-3) 

A i,j 2nr m ‘ Ar 

Referring to the temperature calculation input- instructions, 
card 7, but recalling L must be input in mm not m. 



In the case of heat transfer in the radial direction, 
expression 

H; L, = r m ; L 2 = H 


A- j is obtained from the expression 


A. • « 2Hr 
l,j m 


and similarly for the radiation term abo'»/e 


a *m 




L 2 = 2H 

in which H is the length of the cylindrical surface; where heat 
is conducted between i and j, r m is given by the same equation 
as above (Fig. B-4 (a)); where heat is convected between i and j, 
r m is the radius of the cylindrical surface (Fig. B-4(b)) ; where 
heat is radiated between i and j, t is the radius of the enclosed 
cylindrical surface and r m * the radius of the surrounding cylin¬ 
drical surface (Fig* B-4(c)). 
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B. 2 TRANSIENT ANALYSIS 

For the transient analysis all of the data pertaining to 
the node to node heat transfer coefficients must be provided by 
the input. Additionally, the volume and the specific heat at 
each node is required. For metal nodes this input is straight-¬ 
forward. However, when fluid flow is being considered there is 
no easy way to approximate the fluid nodal volume in a free 
space such as the bearing cavity. However, through use of CYBEAN 
the user's ability to make appropriate estimates will improve. 
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